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Abstract 

We study bounds on Higgs boson masses from perturbative unitarity in the Georgi-Machacek 
model, whose Higgs sector is composed of a scalar isospin doublet, a real and a complex 
isospin triplet fields. This model can be compatible with the electroweak precision data 
without fine tuning because of the imposed global SU (2)r symmetry in the Higgs potential, 
by which the electroweak rho parameter is unity at the tree level. All possible two-body 
elastic-scattering channels are taken into account to evaluate the S-wave amplitude matrix, 
and then the condition of perturbative unitarity is imposed on the eigenvalues to obtain 
constraint on the Higgs parameters. Masses of all scalar bosons turn out to be bounded from 
above, some of which receive more strict upper bounds as compared to that in the standard 
model (712 GeV). In particular, the upper bound of the lightest scalar boson, whatever it 
would be, is about 270 GeV. 
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1 Introduction 



The nature of electroweak symmetry breaking remains unknown at the present status of our knowl- 
edge for high energy physics. In the standard model (SM), a scalar isospin doublet field, the Higgs 
field, is introduced to be responsible for spontaneous breakdown of electroweak gauge symmetry. 
Its vacuum expectation value (VEV) triggers the symmetry breaking, so that it provides origins 
of masses of weak bosons via the Higgs mechanism, and also does those of quarks and charged 
leptons via Yukawa interaction. Although the SM Higgs sector is simple, the Higgs sector could 
have a more complicated structure in the actual world. In particular, when the Higgs sector would 
play an additional role to explain phenomena which the SM cannot, it should necessarily be an 
extended form from the SM one. Therefore, experimental detection of the Higgs boson and pre- 
cision measurements of its properties are extremely important not only to confirm our basic idea 
of electroweak symmetry breaking but also to determine details of the Higgs sector and further to 
outline the structure of new physics. 

In constructing an extended Higgs sector, there are two important requirements from current 
experimental data. First of all, the data indicate that the electroweak rho parameter (p) is very 
close to unity. Second, flavor of quarks and charged leptons is (approximately) conserved in 
the neutral current. In the SM, these two conditions are satisfied respectively by the custodial 
symmetry which ensures p = 1 at the tree level, and by the Glashow-Iliopoulos-Maiani (GIM) 
mechanism which prohibits the tree-level flavor changing neutral current (FCNC). Needless to say 
that these experimental requirements must be respected in extended Higgs models which would 
appear in the low energy effective theory of a more fundamental theory beyond the SM. 

Extension of the SM Higgs sector can be considered by including additional scalar isospin 
singlets, doublets and higher multiplets. It is known that additional singlets and doublets keep 
p = 1 at the tree level [1] . Radiative corrections can slightly deviate the rho parameter from unity, 
corresponding to explicit violation of the custodial symmetry in the dynamics in the loop. On the 
other hand, extension with higher multiplets such as triplets is usually problematic, predicting the 
rho parameter to be explicitly different from unity already at the tree level [2]. One way to avoid 
this problem is to make a fine-tuning on the size of vacuum expectation values of the triplet fields; 
i.e., to set tiny values on them. Another possibility is to impose the custodial symmetry to the 
Higgs sector, so that the rho parameter is predicted to be unity at the tree level. In 1985 Georgi 
and Machacek proposed such a model with one real triplet (Y=0) and one complex triplet (Y=2) 
in addition to the Higgs doublet [3]. Chanowitz and Golden have explicitly constructed the Higgs 
potential of this model [4]; i.e., imposing the custodial SU(2)y symmetry to the potential, VEVs 
of all the isospin triplets become common, and then the tree-level value of the rho parameter 
is unity. They also have shown that the quantum correction from the scalar sector is stabilized 
by such a global symmetry, so that the rho parameter is corrected at the loop level only due 
to explicit SU(2)y violation in the other sectors such as hypercharge interaction and Yukawa 
interaction, just like in the SM. Several phenomenological studies have been done on this model 
in Refs. [5, 6, 7, 8, 9, 10, 11, 12]. 

Generally in extended Higgs models, there are many free parameters in the Higgs potential, 
which spoil predictive power of the model. Hence, it is important to clarify allowed regions in the 
parameter space not only by using current experimental data but also by investigating theoretical 
consistencies such as perturbative unitarity [13, 14], vacuum stability and triviality [15]. This kind 
of study has been often developed to constrain parameters of the Higgs sector in the context of 



the two-Higgs-doublet model [16, 17, 18, 19], and in a specific triplet model [20]. 

In this paper, we study bounds on Higgs boson masses from perturbative unitarity in the 
Georgi-Machacek (GM) model. The Higgs potential respects the global SU(2) R symmetry, so 
that the custodial SU (2) v symmetry remains after the electroweak symmetry breaking (SU(2) L ® 
SU(2)r — > SU(2)y). There are ten physical scalar states, which can be expressed by a SU(2)y 5- 
plet (F 5 ++ , # 5 + , Hi H5, Hs~), a 3-plet (# 3 + , Hi # 3 ~) and two singlets H° and H? [3]. The scalar 
components in the same multiplet are degenerate in mass at the tree level. In the Higgs potential 
of the GM model, explicit Z 2 violation can only appear in the trilinear scalar interaction, but they 
must be forbidden to avoid excessive magnitudes for masses of neutrinos. Neglecting such terms 
by imposing the Z 2 symmetry, all Higgs boson masses in this model are described in terms of the 
VEV, the mixing angles and the dimension-less coupling constants Aj in the Higgs potential. This 
situation is somewhat similar to the two-Higgs-doublet model with the discrete Z 2 symmetry [21], 
in which perturbative unitarity gives upper bounds on all the Higgs boson masses [16]. 

In our analysis, all possible two-body elastic scattering channels (91-channels) are taken into 
account to evaluate the S-wave amplitude matrix in the GM model. Constraints on the Higgs 
parameters are obtained by imposing the condition of partial wave unitarity on the eigenmatrix 
of the S-wave amplitude. Masses of all scalar bosons turn out to be bounded from above, some 
of which receive much stronger bounds as compared to that in the SM (712 GeV). In particular, 
the mass of at least one of the charged Higgs bosons should be less than about 400 GeV. At least 
one of the neutral Higgs boson is lighter than 322 GeV. Furthermore, the upper bound of the 
lightest scalar boson, whatever it would be, can be about 269 GeV. We also find that by using the 
experimental constraints from Zbb results [11], the combined upper bound for the lightest Higgs 
boson is lower than 269 GeV, depending on what the lighest is. Therefore, the model can be well 
testable at current and future collider experiments. 

In Sec. 2, a brief review of the GM model is given. The transition matrix for two-body elastic 
scatterings is calculated in the high-energy limit, and its eigenmatrix is obtained in Sec. 3. In 
Sec. 4, the condition of S-wave unitarity is imposed for the eigenmatrix of the transition matrix, 
and bounds on the Higgs boson masses are evaluated. Conclusions are presented in Sec. 5. 

2 Georgi-Machacek Model 

The GM model contains a complex SU(2) L doublet field (V=l), a real SU(2) L triplet field f 
(Y=0) and a complex SU(2)l triplet field x (X = 2) [3], and respects the global SU(2)r symmetry 
in the Higgs potential [4] . They can be described by the form of SU(2)l <8> SU(2)r multiplets $ 
and A in the potential; 




(2.1) 
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+/Ul Tr($t|$|)A 1 / + // 2 Tr ( A^Tj AT}) , (2.2) 
where r,- are the 2x2 Pauli matrices and 



P f AP, P= 1. (2.3) 




The neutral components of the doublet and the real and the complex triplets have the VEVs, 
vq, V£, and v x , respectively, which are defined as 

0° = ^±^M, (2.4) 

e° = ve + e, (2.5) 

X° = -x + ^^- (2-6) 

After electroweak symmetry breaking, the custodial SU (2)v symmetry remains in the Higgs sector 
(SU(2)l x SU(2)r — > SU(2)y), by which the real and complex triplets have the same VEV, 
va = v$ — v x . Consequently this leads to p = 1 at the tree level [4]. In this case the VEVs are 
constrained as v 2 = v | + 8v A , where v = (y/2G F )~2 ~ 246 GeV. Therefore, differently from usual 
triplet models, i>a can be of order 100 GeV in this model without explicit inconsistency with the 
experimental value of the rho parameter. It is convenient to introduce the doublet-triplet mixing 
angle 9 H , 

tan^ = ^. (2.7) 

Vcj, 

The experimental constraint on 6 H is discussed in Ref. [11]. 

In the potential Eq.(2.2), the last two terms with the coupling constants /ii and /x 2 explicitly 
violate the discrete Z 2 symmetry under the transformation of $ — > $ and A — > —A. Without 
the Z 2 symmetry, the model is allowed to have the mass terms for the neutrinos by assigning of 
lepton number —2 to the complex triplet field, 

i(h v ) ah ^l a CT 2 xi>Lb + ^., (2.8) 

where x — y(-P^x)*- I n order to generate the tiny neutrino masses the Yukawa coupling h u should 
be fine-tuned to be very small as h u ~ O(l0~ 12 ) for the triplet VEV of order 100 GeV. Since we 
would like to avoid such fine tuning with respect to the neutrino masses, we require the discrete Z 2 
symmetry in the Higgs potential and prohibit the last two terms in Eq.(2.2) 1 . Therefore, quarks 
and leptons couple to the SU{2)l doublet field $ in the same way as the SM Yukawa coupling, 
but do not to the triplet A at the tree level. Because all the masses of quarks and leptons are 
obtained from the VEV in we do not have to worry about FCNC, and it is expected to appear 
at most at the same level as that in the SM. This property of the coupling to fermions would give 
an additional constraint on the value of the doublet-triplet mixing angle 8 H by tan 6^ < 0(1), 
since large values of tan^ (3> 1) imply that the top- Yukawa coupling is much greater than 0(1). 



1 The neutrino masses might be generated by any other mechanism (e.g. [22]). We shall discuss it elsewhere [23]. 
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In the GM model, there are ten physical states in the Higgs sector, which are classified as a 
5-plet (H^ + ,H^,HlH^,H^-), a 3-plet (Hf,Hg,Ha), and two singlets H° and H? under the 
custodial SU (2) v symmetry. These are given in terms of the original component fields and the 
doublet-triplet mixing angle Oh as [5] 

(2.9) 
(2.10) 
(2.11) 

W + , (2.12) 

(2.13) 
(2.14) 
(2.15) 

The 5-plet components do not include the component fields from the isospin doublet field $, so 
that the states of the 5-plet do not couple to the fermions at the tree level. On the other hand, 
the 3-plet fields can couple to the fermions. Because of invariance under the custodial SU(2)y 
symmetry, states in the different multiplet cannot mix each other. 

All members in the same SU(2)y multiplet are degenerate in mass at the tree level. The 
masses of the 5-plet and the 3-plet are respectively given by 





+ = X ++ , 




Ht 


= ( X + -n/v / 2, 




Hi 


= (2C° - V2 X ° r )/VQ , 




Ht 


= cos9 H ( X + + e)/V2- 


— sin 


Hi 


= i(- cos HXi + sin9 H c 




Hi 






H°i 


= (v / 2x° + C°)/v / 3. 





m 2 H5 = (X 4 sm 2 9 H - ^X 5 cos 2 H )v 2 , (2.16) 



< = -y- 2 - (2-17) 



On the other hand, two SU (2)v singlets can mix, and the mass matrix 



Scos 2 ^ 4sin^cosM2A3 + A 5 ) 

1 v /|sin^cos^(2A3 + A 5 ) sin 2 9 H (3X 2 + A 4 ) 



is diagonalized by introducing the mixing angle a. The eigenvalues correspond to the masses m H o 
and m H o> for the mass eigenstates Hf and H± . 

From Eqs.(2.16) - (2.18), the quartic couplings A« are expressed in terms of the masses and the 
mixing angles as 



Ai = (m^ocos a + m H o> sin a)/(8v cos 9 H ) , (2-19) 

A2 = ( m i5-o sin 2 a + rn 2 H o> cos 2 a — m 2 Hs + Sm 2 H3 cos 2 9h) / (3f 2 sin 2 Oh) , (2.20) 

A 3 = (m 2 H o' — m 2 H o) cosa sino;/(v / 6t> 2 sin H cos ^) + rn 2 H3 /v 2 , (2.21) 

A 4 = (m| 5 - 3m^ 3 cos 2 ^)/(^ 2 sin 2 H ) , (2.22) 



A 5 = -2m 2 jv 2 . (2.23) 



The SU (2) v 3-plet fields receive the constraints from the current data of Z — > fefe, _B — -Bo and 
-Ko — -^0 mixings [24, 25]. These data give bounds on the mass m H3 with the mixing angle Oh- 
The most stringent experimental constraint comes from Z — > bb. The mass in Hz is constrained to 
be smaller than 1 (0.5) TeV for tan## < 2 (1). 
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Although the 5-plet fields do not couple to the fermions, the singly-charged state in the 5-plet 
has a characteristic coupling of H^W T Z, which only appears beyond the tree level in multi- 
Higgs-doublet models [26]. Experimental confirmation of a sizable coupling of HfW T Z with 
p ~ 1 should be a strong indication for the GM model [27]. This coupling is testable via the 
process pp — > W ± H T at the Fermilab Tevatron [12], also via pp — > W ± *Z*X — > H ± X [28] and 
the decay process — > W T Z [29] at the CERN LHC, and further via the processes e + e~ — > 
W^if* [8, 9, 30] and e+e~ -> vvW^Z* -> z/z/iy* [31] at the ILC. Another striking feature of 
models with complex isospin-triplets, such as the GM model, the left-right symmetric model, the 
littlest Higgs model, and some models motivated by neutrino masses, is the appearance of doubly- 
charged states i/ ±± . At hadron colliders, such doubly-charged Higgs bosons are studied via the 
pair production mechanism [10, 32, 33] as well as the single production mechanism [33, 34] and 
the lU-boson fusion mechanism [6, 35]. They can also be investigated at the ILC and its e~e~, 
e~7 and 77 option in various scenarios [36]. 

3 The S-matrix for two-body elastic scatterings 

In this section, we calculate the transition matrix of elastic scatterings of two scalar-boson states 
in the GM model. The transition matrix T(tpitf2 — > V 9 3 ( / 9 4) is equivalent to the S-wave amplitude 
( < ^3 < ^4l a °lv 9 iV 9 2) at high energies (\/s ^> m^y), where tpi represent longitudinally-polarized weak 
bosons or physical Higgs bosons of the model. The condition of partial wave unitarity is given for 
the S-wave amplitude matrix by [2, 13] 



We employ this condition in the high energy limit to constrain the model parameters in the next 
section. Thanks to the equivalence theorem [37], the S-matrix elements in which longitudinally- 
polarized weak bosons are in initial and final states are equivalent to those in which these weak 
bosons are replaced by the corresponding Nambu-Goldstone bosons in the high energy limit [13] . In 
addition, in this limit, only quartic couplings (scalar contact interactions) of the Higgs- Goldstone 
couplings are relevant to the unitarity conditions, which can be translated into the bounds on the 
related Higgs-boson masses after Eq.(3.24) is imposed. Therefore, we here evaluate the matrix 
((/93(/9 4 |a°|(/9 1 (/9 2 ) taking into account all possible two-body scalar channels in the high energy limit, 
and obtain all the eigenstates and the eigenvalues. 

Under 0(4) (~ SU (2) L <g> SU(2) R ), the field components of 0, x an d £ are expressed by a 4 and 
a 9 representations as 



where 0+ = (^ + iu 2 )/y/2, 0° = (<f>° r + #°)/\/2, X ++ = (Xi + %)/v^, X + = (Xs + iX*)/V2, 



X° = (Xr +^X?) / and £ + = (£1+^2) / V%- We consider all possible two-body scattering channels 
(^a^b ^c^d) not only for the neutral two-body states as initial and final states but also for the 
singly-, the doubly-, the triply- and the quadruply-charged two-body states. There are totally 91 
initial (or final) two-body states, in which 25 are the neutral, 36 are singly-charged, 22 are doubly- 
charged, 6 are triply-charged, and the last 2 are the quadruply-charged states. We construct the 
91x91 transition matrix of high-energy S-wave amplitudes, and then calculate their eigenvalues. 



|(y?3V?4|a \<pi(p 2 )\ < ~- 



(3.24) 




(3.25) 
(3.26) 
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The initial (final) two-body states can be treated separately as ^d^d, ^t^t and ^d^t- The 
high-energy S-wave amplitudes are block-diagonalized by the electric charge and also the discrete 
Z 2 symmetry ($ — > $ and A — > —A). Each submatrix with respect to the \I/z)\I/_d or ^f T ^f T states 
can also be classified by irreducible decomposition of direct products of the representations for 
0(4) as 

4®4 = (1) D ©(9)©(6), (3.27) 
9® 9 = (l) T © (44) © (36). (3.28) 

The only singlet and symmetric representations, 

= E(^) 2 , (3-29) 

k=l 

(s) ij = *i*i--^E(*2) 2 . ( 3 - 3 °) 
iVa k=i 

contribute to the scatterings of our interests, where a = D or T; i,j = 1—4 and s = 9 for a = 
D , or i, j = 1-9 and s = 44 for a = T, and N D = 4 and iV T = 9. For the ^f D ^f T states which 
should be of the 4® 9 representation, there is no singlet representation so that 0(4) cannot help for 
the classification. Furthermore, several additional discrete transformations can be used to further 
classify the states, which will be defined below. 

Neutral channels 

We outline further classification of the decomposed irreducible states for the case of the neutral 
25 two-body channels (4 for ty D ty D , 11 for ty T ty T and 10 for ty D ty T ). For ^ D ^I D states, we have 
the singlet state (1)^ and the three neutral elements of (9)* J ((9) 33 , (9) 44 and (9) 34 ), in which the C 
parity separates (9) 34 from the other states. After taking appropriate linear combination, we obtain 
two separate states under the transformation of 4> r — > 0; and fa — > —fa as ((9) 33 ± (ty 44 )/^. 
Thus four eigenstates of the transition matrix for the neutral ^d^d — * ^d^d channels are 
obtained [13]. 

Next, we consider ^t^t — ► ^t^t scatterings in which both the initial and final states are 
electrically neutral. In addition to the singlet state (1) T , we have 10 neutral states from (44)* J , in 
which { (44) n +(44) 22 , (44) 33 +(44) 44 , (44) 55 +(44) 66 , (44) 77 , (44) 88 } are the diagonal element states 
(i = j), and {(44) 35 + (44) 46 , (44) 36 - (44) 45 , (44) 78 , (44) 79 , (44) 89 } are the off-diagonal element 
states (i ^ j). Among the diagonal element states, the linear combination (44) 77 — (44) 88 has 
different property under the transformation of Xr — * Xi and Xi ~~ ~Xr- Then, linear combinations 
{ (44) 11 + (44) 22 + (44) 77 + (44) 88 , (44) 33 + (44) 44 + (44) 55 + (44) 66 } and { (44) 11 + (44) 22 - 
(44)77 _ (44)88 ? ^44^33 + (44)44 _ (44)55 _ (44)66 j show different property under the transformation 

of x ++ X ^ and X + X~ ^ i + i ~ ■ The first two states have completely the same property 

as that of the singlet state (1)t, so that the appropriate linear combination of these three states 
give the three eigenstates. For the off-diagonal element states, we can separate them by using the 
C parity and the transformation of £ — > — £, so that these states are block-diagonalized to two 
2x2 submatrices and one singlet. By diagonalizing remained 2x2 matrices, we obtain all the 
eigenstates for the ^t^t — > ^J T ^ T channels. 

In order to diagonalize all the ^d^d and ^t^t states, we take linear combinations of the 
eigenstates of ^d^d and ^/t^t that have similar transformation properties. Consequently, all the 
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typ^D and ^t^t states are completely separated, and the eigenvalues of the transition matrix 
for these channels are obtained. 

Finally, classifying the ten \l/ D ^ T states by using the C transformation as well as some discrete 
transformations in a similar way to above, we completely diagonalized 25 x 25 scattering matrix 
for the electrically neutral two-body states. 

The neutral states, \A®)(i = 1 ~ 25), which give (at most 2x2) block-diagonal transition 
matrices are found as below: 



l^l> 




(20+0- + 0°0° + 0°0°)/2 v / 2 , 


(3.31) 


1 Z / 





(2x ++ X~ + 2x + X~ + 2^ + r + XiXi + X°X° + CC)I^ , 


(3.32) 


\A° 3 ) 





(20+0- - 0°0° - 0°0°)/2v / 2 , 


(3.33) 







(-2 X ++ X~~ + XiXi + XrXr)/lV2 , 


(3.34) 




— 


(0-0- - , 

\ 1 1 1 1, 1 1 1 1 / 1 ' 


(3.35) 


\<) 


= 


(x + c+x-e+V2xX)/2, 


(3.36) 


\A° 7 ) 


= 




(3.37) 


\A°\ 

1^8/ 


= 


(x + r - x-e - v / 2ic o x°)/(20 , 


(3.38) 


\A°\ 


= 


(X?X? - X?X?)/2 , 


(3.39) 


4 U \ 


— 


(2x ++ X~ _ - X + X _ - + X?X? + X?X° - 2C°C)/3V2 , 


(3.40) 






\ A X X 4 X X 4| S S + XrXr + XiXi + 4 Sr-?rJ/ D vZ , 




14° \ 







(3.42) 


14° \ 

1-^13/ 


= 


(xV + x-r-v^^/a, 


(3.43) 


14° \ 


= 


(x + r - X^ + + v / 2^°X°)/(20 , 


(3.44) 


1 A \ 

^15/ 


= 


n n 

V V 
ArAi 5 


(3.45) 


4 U \ 
1 ^16/ 


= 


|2(0 + x" + <Tx + + 0V + <T£ + ) - v/2(0°X° + <t>iXi) ~ 20°C°} /3 , 


(3.46) 


14° \ 
1-^17/ 




(<« - fe )/^ , 


(3.47) 


14° \ 




(0^ + 0°X°)/v / 2, 


(3.48) 


14° ) 




{2(0 + x" - <Tx + - 4> + C + <T£ + ) + v/2 ? (0°x° - 0°X°) - 4^C°} /(«) , 


(3.49) 


14° \ 

1^20/ 




{0 + X~ + <Tx + + 0V + + v / 2(0 r °x r ° + 0°x°) + r °e°} /3 , 


(3.50) 


14° \ 




{0 + X" - 0^X + - + r + + v / 2^(0M ) - 0M + i^r} /(3i) , 


(3.51) 


14° \ 




{<p + x~ - ^x + - + r + + 2^(0^ - 0°x° ) + 4*0°e } /(^ , 


(3.52) 


14° \ 

l /1 23/ 




{0 + X~ + 0"X + + + r + - 2v / 2(0 r °x r ° + 0?x-) + 40 r °C°)} /6 , 


(3.53) 


14° \ 
l /1 24/ 




(0+x- + 0~x + - </> + r - <rn/2 , 


(3.54) 


14° \ 
l /1 25/ 




(0 + x~-0~x + + + r-0^ + )/(2O . 


(3.55) 



The state (l)^ and the three neutral linear-combined states from (9) iJ respectively correspond to 
| A?) and {\A° 3 }, \A° 5 ) and \A%)}. The state (1) T and the linear combined states from (U) ij corre- 
spond to \A®) and \A®) (i = 4,6,8, 11 — 15), respectively. The first eight states block-diagonalize 
the transition matrix to four 2x2 submatrices, and the other seventeen states give eigenstates. 
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Singly- charged channels 



There are eighteen singly-charged states with the electric charge +1 among all the two-body 
states. The charge conservation ensures that these states are composed of a subset among all the 
states with various electric charges. The corresponding high-energy transition matrix is conse- 
quently (block-) diagonalized by the following states: 



1 1 / 


= (0 + 0° + + 0°)/v / 2, 


(3.56) 


\M) 


= (V2 X ++ X~ + X + Xr + X + X?) /2 , 


(3.57) 


\M) 

1 o / 


= (0+0° - 0+0°)/^ , 


(3.58) 


K> 

i *± i 


= (V2 X ++ r + ex°r - e + x°)/2 , 


(3.59) 


1 / 


= (V2 X ++ X~ ~ X + X° r - X + X°i + 2£ + C°)/2\/2 , 


(3.60) 


1 / 


= (v / 2x ++ r - t + Xr + £ + X? + 2x + £°) /2v^ , 


(3.61) 


\M) 

\ ( I 


= (v / 2x ++ X _ - X + X? - X + X? - 2£ + C°) /2\/2 , 


(3.62) 


\At) 

1 o / 


= (v^x^r - £ + Xr + £ + X? - 2x + C°)/2v / 2 , 


(3.63) 


\Ai) 

\ y / 


= (£ + X? + £ + X?)/v / 2 , 


(3.64) 




= (x + x r ° - X + X°)/V2 , 


(3.65) 


I4i> 


= (x ++ 0~ - 20+C° " X + 0° - X + 0° - t€ + e + 0°)/3 , 


(3.66) 


\Af 2 ) 


= (0 + X° + + X° + v / 2e + 0° + v / 2e + 0°)/v / 6, 


(3.67) 




= (0 + X° - + X° + v/2x + 0° - v / 2x + 0°)/v / 6 , 


(3.68) 




= (4 X ++ 0- + 40+ Xr ° - X +0 r ° - X + 0° - t€ + t<P- )/6 , 


(3.69) 




= (x + 0° + X + 0°"« + e + 0°)/2, 


(3.70) 


\At 6 ) 


= (v^+x? + - « - f^/Ve , 


(3.71) 


\Atr) 


= (v / 20 + X r °-V20 + x°-X + ^° + x + 0-)/V6, 


(3.72) 


K 8 > 


= (2 X ++ 0- - + C° + X + 0° + X + 0? + - e + 0°)/3 • 


(3.73) 



The eighteen singly-charged states with the electric charge —1 can be obtained by the C trans- 
formation for the above states with the charge +1. 

Doubly- charged channels 

There are eleven doubly-charged two-body states with the electric charge +2. We can decom- 
pose the subset of the transition matrix for these states to at most 2x2 matrices by the following 
linear combination; 



\A+ + ) = 0+0 + , (3.74) 

\At + ) = (x + e-X ++ e r )/V2, (3.75) 

l^ 3 ++ ) = X + X + , (3.76) 

\At + ) = (X ++ X° r -X ++ Xf)/V2, (3.77) 

\Aj + ) = (3.78) 

\At + ) = (X ++ X°r + X ++ X°)/V2, (3.79) 

\A+ + ) = ( X + e + X ++ £)/V2, (3.80) 
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K + > = (0 + X + + + e + + X ++ 0°)/v/3, (3.81) 

K + > = (0 + X + -0 + £ + -X ++ $)/V3, (3.82) 

Ko + > = (0 + X + + 4> + t ~ 2x ++ €)/V6 , (3.83) 

Ki + > = (0 + X + -^ + + 2x ++ 0°)/v / 6. (3.84) 

The corresponding doubly-charged two-body states with the charge —2 can be obtained by C 
transformation in the above states with the charge +2. 

Triply-charged channels 

There are three triply-charged states with the electric charge +3, and the subset of the tran- 
sition matrix for the initial and final states can be diagonalized by the following eigenstates as 

\A+ ++ ) = (3.85) 
|^ ++ ) = (3.86) 
\A+ ++ ) = X ++ e- (3.87) 

All the eigenstates with the opposite electric charge can be obtained by the C transformation of 
these eigenstates with the charge +3. 

Quadruply-charged channels 

Finally, we have only one quadruply-charged state for each electric charge of +4 and —4, 

]At+++) = (3 _ 8g) 

l^r ) = X~~X~~- (3-89) 
Eigenvalues of the transition matrix for all channels 



In summary, the transition matrix T has been block-diagonalized as 
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(3.90) 



where block-diagonal transition submatrices for the neutral, the singly- charged, the doubly-charged, 
the triply-charged, and the quadruply-charged two-body states, T°, T 1 * 1 , T ±± , T ±±=t , and T ±±±=t , 
respectively, are given by 

T° = diag(Xi, X 2 , X 2 , X 2 , y u y u y 2 , y 2 , 2/2, 2/2, 2/2, 2/3, 2/3, 2/3, 2/3, 2/3, 2/4, 2/4, 2/4, 2/5, 2/5) , (3.91) 
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= diag(X3,X3,y 6 ,y 6 ,y 7 ,y 7 ,y 2 ,y 2 ,2/3 ) 2/3,2/3,2/4,2/4,2/4,2/4,2/5) , (3.92) 
T ±=t = diag(X4,X5,X 5 ,?/2,2/3,2/3,2/4,2/4) , (3.93) 
T ±±± = diagd/3,2/2,2/2) , (3.94) 



Here Xi are the 2x2 matrices whose eigenvalues xf are given by 



= 12Ai + 22A 2 + 14A 4 ± v^ 12 ^! ~ 22A 2 - 14A 4 ) 2 + 144A§ , (3.96) 
= 4Ai + 4A 2 - 2A 4 ± V / (4A 1 - 4A 2 + 2A 4 ) 2 + 4A§ , (3.97) 
= 4A 2 + 4Ai ± V /(4A 2 - 4A X ) 2 + 4A| , (3.98) 



X 2 



J- n 



xf = 8Ai+4A 2 -2A 4 ± V(8Ai -4A 2 + 2A 4 ) 2 + 8Ai , (3.99) 



xf = 12A 2 + 14A 4 ± 2 y/4Ai + 4A 2 A 4 + 17A| . (3.100) 
The eigenvalues y^ are obtained as 

yi = 8A 2 + 16A 4 , (3.101) 

y 2 = 8A 2 + 4A 4 , (3.102) 

y 3 = 4A 3 + A 5 , (3.103) 

y A = 4A 3 -2A 5 , (3.104) 

y 5 = 4(A 3 + A 5 ) , (3.105) 

y 6 = 8A 2 + 4(2 + V2)A 4 , (3.106) 

y 7 = 8A 2 + 4(2 - V / 2)A 4 . (3.107) 

Although the transition matrix between initial and final two-body states is originally 91 x 91, the 
number of independent eigenvalues turns out to be only seventeen. 

4 Unitarity bounds on the masses 

In this section we analyze mass bounds on the Higgs bosons in the GM model, imposing the condi- 
tion of perturbative unitarity in Eq. (3.24) to the transition matrix given in the previous section. 
Consequently we obtain seventeen inequations with respect to all the independent eigenvalues of 
the transition matrix T in Eq.(3.90) as 

\ x tl \ x f\, \ x i\, \ x 4\, \ x t\, \yi\, I2/2I , I2/3I, M, |a/s I , |ye|, \y7\ < 8vr. (4.108) 

These eigenvalues are respectively given in Eqs. (3.96) - (3.107) as a combination of the dimen- 
sionless coupling constants Aj (i — 1 ~ 5) in the Higgs potential, and Aj are related to the Higgs 
boson masses through Eqs.(2.19)-(2.23), these constraints can be translated into the bounds on 
the masses rrijjo, mg i ,m H3 and m Hs and on the mixing angles 6 H and a. 

We here show the numerical results on the Higgs mass bounds. Fig. 1 shows the allowed 
regions of the masses in the rn H3 - rn H5 plane (a), in the m H3 - m^o plane (b) and in the m^ 
- rrijjo' plane (c). We vary the Higgs boson masses in the range % 3 ,% 5 ,%o,%o' < 1 TeV 
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Figure 1: Allowed regions of the masses of the Higgs bosons in the m H:i - m H5 plane (a), in the 
m H 3 ' m H° plane (b) and in the m H o - m H o> plane (c). In Fig. 1(a) the light shadowed regions are 
excluded by the Zbb results. 



and the mixing angles for < Oh < vr/2 and — n/2 < a < 11/ 2. In each figure, the conditions of 
perturbative unitarity in Eq.(3.24) are satisfied inside the regions. In Fig. 1(a), light shadowed 
region is excluded by the Z — > bb result. 

Fig. 1(a) shows that m H and m Hr are bounded from above respectively by about 400 GeV 
and about 700 GeV. These upper bound come from \xf\ < 8n, which give the most stringent 
constraint among the inequations in Eq.(4.108). For m H3 ^ 170 GeV, is bounded from above 
whose border is approximately corresponding to A 4 ~ or m Hh ~ \f%m Hz with 9 H ~ 0. Owing to 
the factor \/3, the more strict constraint is given on m Ha than on m H . For m Hz <^ 170 GeV, on the 
other hand, the upper bound on m Hz> (290 GeV) is realized at 9 H ~ ir/2. When all masses other 

than m Hz are zero, m H3 is bounded by m,H 3 < \J2tt/3v (~ 356 GeV) from \xf\ < 8n. However 
numerical analysis shows that the actual upper bound is a few decade GeV greater. This excess 
comes from some delicate cancellation in Eq.(3.96). When we impose the experimental data from 
Z — > bb which give the constraint on the combination of m H3 and 9 H [11], the allowed region is 
further limited in the dark shadowed regions. The upper bounds on m H3 and m H5 do not change, 
but the remained allowed region is in the vicinity of rn H5 ~ V3m H3 . 

In Fig. 1(b) we can see that the upper bound on m H o is about 710 GeV, which is almost 
the same as that on the mass of the SM Higgs boson [13]. Larger values of m H o are allowed 
for smaller values of a and 9 as well as smaller m H5 and ttl h o' values. For instance, taking the 

limit a — > 0, 9 H — > 0, m Hr> — > and m H o> — > 0, we obtain m H o < ^8ii/3v (~ 712 GeV) in 

the condition \xf\ < 8n. The mass bound for another singlet Hf' can be obtained by replacing 
m H o with m H o> and a with a + n/2, which can be seen from Eqs.(2.19)-(2.21). Consequently 
the allowed regions in the rn Hz - m H o> plane are given by the same as in Fig. 1(b). We find that 
contrary to the result in Fig. 1(a) there is only few difference in the case where we include the 
Z — > bb data. 

The allowed region in the m H o - m H o> plane in Fig. 1(c) is symmetrical about the line of 
m H o = m H o'(= m). It is notable that at least one singlet receives very strict constraint from 
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perturbative unitarity. The mass of lighter singlet, either H® or H®' , is bounded from above by 
322 GeV. In analytic calculation, this upper bound is obtained as m < y/6ir/llv (~ 322 GeV) 
from \xf | < 8tt. 

In the following, we evaluate the upper bound on m lightcst = Min (jn Hs , m H5 , m H o, n^jjo'^j ■ We 
start from the case in which the constraints from the Zbb results are switched off. When the 3-plet 
is the lightest, we obtain the upper bound on mutest (= m H 3 ) as 

^lightest < 269 GeV, (4.109) 

which is considerably lower than that of the SM Higgs boson, 712 GeV. This condition comes from 
the constraints \xf\ < 8n and \x$\ < 8ir. Similar analysis has been done for m lightcst = m Hs ,m H o 
and rrifjoi in order, and the same bound as in Eq.(4.109) is derived for each case. When m u htest ~ 
269 GeV, all the masses are degenerate in mass {m H:j = vn Hh — m H o = m H o'). The situation 
turns out to be quite similar to the situation of the two-Higgs-doublet model with the discrete 
symmetry, where the lightest of all Higgs masses are bounded at 410 GeV [16]. In the case of the 
GM model the number of the two-body states is greater than that in the two-Higgs-doublet model. 
(The neutral two body states are 14 channels in the two-Higgs-doublet model and 25 channels 
in the GM model.) Thereby we have obtained the stronger bounds than the two-Higgs-doublet 
model. Finally, when we take into account the Zbb results [11], the angle 6h is more limited for 
smaller values of rn Hr Consequently, the combined upper bound on m lightest becomes lower than 
269 GeV. Depending on what the lightest is, the combined upper bound turns out to be about 
249 GeV (176 GeV) when m Hs , m H o or m H o> (m Hs ) is the lightest. 

We have not included the LEP direct search results, which give the lower bound itltt > 114 
GeV in the SM [38]. In the GM model similar lower mass bounds can be obtained for neutral 
Higgs bosons but depending on the mixing angles, which would slightly affect the upper bounds 
by using the results in Figs. 1(a), (b) and (c). We have taken into account only the Z — > bb result 
as the experimental constraint [11], because this constraint drastically changes the bound in the 
m H s ~ m H 5 plane and also that on mutest- 



5 Conclusions 

In this paper, we have analyzed unitarity constraints on the Higgs boson masses in the GM 
model, which includes a real and a complex isospin triplet fields but predicts p = 1 at tree level. 
All possible two-body elastic-scattering channels (91 channels) have been taken into account to 
construct the S-wave amplitude matrix in the high energy limit. The condition of S-wave unitarity 
in Eq.(4.108) has been applied to the eigenmatrix. 

We have found that all the Higgs bosons receive their masses from the VEV under the discrete 
Z2 symmetry, so that all the masses can be bounded from above by the condition of parturbative 
unitarity. In particular, the upper bound on the mass of the SU(2)y 3-plet is about l/\/3 lower 
than that on the SM Higgs boson mass (712 GeV). Hence at least one of the singly- charged 
Higgs boson masses is bounded from above at about 400 GeV. The mass of the lighter SU(2)y 
singlet scalar state, either Hq or H' Q) turns out to be bounded from above by about 300 GeV. 
Furthermore, the mass of the lightest Higgs boson among the 5-plet, the 3-plet and the two singlet 
states, whatever it would be, receives very strong constraint from above; i.e., Tightest — 270 GeV. 
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The point of the parameter space at which m u ghtest takes its maximum value corresponds to that 
where all the mass parameters are degenerate. The combined upper bound with the Zbb results 
becomes about 150 GeV (95% C.L.). Therefore, the model turns out to be well testable at collider 
experiments. The SU(2)y 5-plet and 3-plet have the doubly- and singly-charged states, so that 
the distinctive phenomenological features of this model should also appear in physics of charged 
Higgs bosons. Detailed phenomenological features will be discussed elsewhere. 

In the analysis above, we have considered the Higgs potential with the Z 2 symmetry, neglecting 
the trilinear scalar terms of /ii and /i 2 . The imposition of the Z 2 symmetry in our analysis would 
be justified to avoid large excess of the neutrino masses. When we do not respect the Z 2 symmetry, 
the upper bounds in above results become relaxed according to the scales of Hi and fi 2 which have 
linear mass dimension. Unless fi± and fi 2 are substantially larger than (9(100) GeV, our results 
above can sufficiently be applied by small relaxation. 

Finally, in this paper, we have employed partial wave unitarity to constrain parameters of the 
GM model at the tree level. A more detailed study with the radiative effects such as vacuum 
stability or triviality might give more strict bounds on the Higgs boson masses in this model. 
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